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three-level A medium, in contrast to previous treastments that have considered the cases of media whose 
ground states are characterized by probabilities (level populations) or by probability amplitudes (coherent 



lO ' pure states). We present analytic solutions of the Maxwell-Bloch equations, and we extend our analysis 

(N : 

with numerical solutions to the same equations. We interpret these solutions in the bright/dark dressed 
r^ ■ state basis, and show that they describe a population transfer between the bright and dark state. For 

mixed-state A media with partial ground state phase coherence the dark state can never be fully populated. 
5 ' This has implications for phase-coherent effects such as pulse matching, coherent population trapping, 

3 , and electromagnetically induced transparency (EIT). We show that for partially phase-coherent three- level 

media, self induced transparency (SIT) dominates EIT and our results suggest a corresponding three-level 



area theorem. 



INTRODUCTION 



The description of radiative phenomena in terms of intensities and probabihties is often satisfac- 
tory but the effects of wave coherences are then neglected. When coherence effects are prominent 
a wave description is usually adopted. Propagation of short laser pulses in resonant media often 
falls into an intermediate domain where neither description is satisfactory. The most important 
time scale is much longer than the period of the laser field but shorter than the decoherence time 
of the medium, and the slowly varying envelope and rotating wave approximations are engaged to 
simplify the description of evolution in this domain [l|, y, y] . 

However, theories of resonant propagation usually ignore the very real possibility that the state 
of the medium may be prepared in a way not adequately described by probability (population) 
assignments to the levels. The term "phaseonium" was introduced by Scully [4| to describe a three- 
level atomic medium in the A configuration, where two ground levels of the atoms are prepared 



in a phase-coherent pure-state superposition. Phase coherent effects lead to coherent 
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the d«k »ea .h»rem U . simul.o.. [14], and adiabatons mM- These effects a.e a. ah governed 

by dark-state [17] considerations. 

Pure phaseonium is not easily prepared experimentally, and the question what effects may 
acompany more realistic preparation has remained open. Here we report new results obtained 
for propagation in media that could be called partial phaseonium, or "mixonium", because the 
medium is allowed to have mixed state or partial coherence in its ground state. The simple cases 
we study are distinct from traditional EIT scenarios in that the ground state coherence is prepared 
ahead of any pulses entering the medium, while in traditional EIT it is the pump and probe pulses 
which prepare the medium coherence. Methods to prepare a phaseonium medium have been shown 
for ladder UM and A 191 systems. 

fT 

A previous study [2u\ by Kozlov and Eberly of two-pulse propagation in absorbing media made 
use of the Maxwell-Bloch model (see e.g. ll|), which describes fully coherent pulse propagation 
through atomic media. They reported two types of pulse evolution, one similar to that of one- 



pulse self-induced transparency (SIT) [2l|], and the other showing a different form of two-pulse 
evolution under conditions they identified with EIT, in which the dark state plays an important 
role. They defined the SIT-type propagation in three-level media as occurring when the dark 
state was nearly or completely empty and EIT-type propagation as occurring when the dark state 
was highly populated. They showed that EIT-type dominated SIT-type for propagation in a fully 



coherent phaseonium medium. However, when a A medium is prepared without ground-state 
phase coherence, the role of the dark state is minimized, and we have shown 23,123] that SIT-type 
propagation is dominant. 

In reahstic experimental preparation, pure-state three-level phaseonium is difficult to achieve. 
Ground state phase coherence can be lost due to environmental decoherence or imperfect prepara- 
tion techniques. Thus we report new results on two-pulse propagation through a A medium that 
is prepared with the ground states in a partially phase-coherent superposition, between the two 
extremes of completely mixed or completely pure state, the basis for our term "mixonium" . We 
include a variable parameter A in our initial state definition that will allow us to examine the prop- 
agation dynamics both analytically and numerically over the entire range from completely mixed 
to completely pure. We will follow the previous definition and determine whether it is EIT or SIT 
that dominates pulse propagation in mixonium. 

We will use a three-regime language to distinguish three propagation zones, similar to our 
analysis of two-pulse propagation in a completely mixed-state medium [23] . In the pure phaseonium 
case the analytic solutions to the Maxwel-Bloch (MB) equations predict that simulton pulses 



14l ] that begin entirely in the strongly interacting bright state will be transferred into the non- 
interacting dark state just as one would expect based on coherent population trapping and dark 
area theorem arguments. If the medium is initially prepared only partially phase-coherent, a similar 
simulton transfer process occurs. However the lack of complete phase coherence prohibits the dark 
state from being fully populated, and the pulses will continue to interact with the ground to excited 
state transition. Our analytic solutions allow us to determine the maximum population of the dark 
state, when the medium is not in a pure state. 

The inability of the dark state to be fully populated, and the subsequent pulse-medium interac- 
tion that continues for non-pure medium preparation causes SIT-like effects and reduces the role of 
EIT. This has dramatic consequences for the pulse matching that occurs in phaseonium. We find 
that our analytic solutions have broad predictive capabilities for mixonium due to pulse-resha ping 
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caused by SIT. In the phaseonium case our numerical solutions agree with previous studies 
predicting temporal pulse matching. However, in the mixonium case SIT begins to play a role 
because the dark state cannot be fully populated. We find that for long propagation distances SIT 
always dominates in a mixed-state medium. Further highlighting this effect, we also show that the 
bright pulse area behaves in a manner very similar to SIT, which leads us to believe that there is 
an equivalent three-level area theorem. 



II. PHYSICAL MODEL 




FIG. 1: Three level atom in the A configuration with level 1 connected to level 3 via a laser field fia, which 
we refer to as the pump field, and level 2 connected to level 3 with laser field fib which we refer to as the 
Stokes field. We assume a two photon resonance condition so that both fields are detuned from resonance 
by an equal amount A. Loss from the excited state is included by a damping term 73. 



We consider dual-pulse propagation in a medium of three-level atoms in the lambda configu- 
ration as shown in Fig. [H Possible physical realizations appear to be in a D line of cesium or 
rubidium, and we give some of the experimental parameters at the end of this section. We assume 
a two-photon resonance condition such that both laser fields are detuned by an equal amount A. 
The Hamiltonian of the system in the rotating wave picture is given by 



H = hA\3){3\ - n— 11)(3| - h—\2){3\ - h^\3){l\ - h^\3){2\, 



(1) 



where ila = 2(ii -Sa/^ is the Rabi frequency of the pump laser field. Throughout this paper, we will 
refer to Qa as the pump field and il^ as the Stokes field to be consistent with the common nomen- 
clature associated with stimulated Raman scattering. Our notation implies that the laser fields can 
be represented as a slowly varying envelope times a carrier frequency, e.g., Ea = £^e'^(^ax-ujat) _[_(;. g_^ 
where ka and uJa are the wavenumber and carrier frequency of the pulse and £a is the envelope 
function. The term di is the dipole moment of the 1 — > 3 transition. Similar notation applies to 
the Stokes Rabi frequency $7;,. 

Individual density matrix equations can be derived from the Hamiltonian and the von Neumann 



equation ihdp/dT = [H, p] . They are given by: 

-Q^ = «^P3i - «^Pi3 (2a) 

-Qf- = « Y^32 - I^P23 (2b) 

-^jT = -«^P31 + ^^Pl3 - « Y^32 + «^/323 (2c) 

dpu .^a .^l ,„,^ 

-^ = «^/>32 - « Y^^^ ^ ^ 

-^ = ZA/913 - ^^Pl2 + «^(/533 - Pll) (2e) 

-^ = zA/923 - «^/'21 + 1^{P33 - P22)- (2f) 

We assume that the temporal duration of each laser pulse envelope is sufficiently short that we 
can neglect decay terms, such as arising from loss to other atomic states, spontaneous emission, 
or coUisional dephasing effects. For alkali vapors this requires pulse durations on the order of, or 
shorter than, 1 ns. By making use of the slowly varying envelope approximation and the rotating 
wave approximation we can reduce Maxwell's wave equation to two independent first order wave 
equations for each individual pulse. They are given by: 

dQa 



dZ 



/oo 
dAF{A)pi3 = -ipaipis) (3a) 

'OO 

/oo 
dAF{A)p23 = -ipb{p23), (3b) 

'OO 



dZ 

where we have written equations ^ and ([3]) in a retarded time coordinate system such that T = 
t — x/c and Z = x/c. Thus the derivatives are given by: d/dt = d/dT and d/dZ = cd/dx + d/dt. 
We use the bracket notation to symbolize a statistical average to take into account inhomo- 
geneous broadening, for example due to thermal motion of the atoms if the medium is a vapor. 
When needed, the average is performed with the function F{A) = T^/v^Tre"' 2) (^-^) /2^ where 
T2 is the inhomogeneous lifetime and A is the detuning of the laser fields from atomic line center. 
We will consistently assume line-center tuning, so A = is implied throughout. The parameters 
Pa = Nd\ijJa/heQ and p^ = A^d^w^/^ieo, where A^ is the density of the atoms, are proportional to 
the usual attenuation coefficient or inverse Beer's length: 



QD(A) = 7rF(A)///c ^ ai5(0) = V^TaV/c, (4) 

for each transition, where the simplified final form applies to line-center tuning. When pa = Pb = P, 
which we will assume hereafter, Eqs. ([2]) and ([3]) are exactly solvable by using methods such as 



inverse scattering 



24 



25l | or Backlund transformations 26|, |27|, [28 1 



A single weak pulse presents a familiar case, in which the pulse causes some population ex- 
change between ground and excited states, and even though we neglect homogeneous loss terms, 
Doppler (or any other inhomogeneous) broadening serves as a dephasing mechanism. Thus a single 
weak pulse will be absorbed as it propagates, and its peak intensity will decay exponentially as 
|r2(x)p = |il(0)pe~'^-^^, where Q,{0) is the peak intensity of the weak input pulse and a^ is the 
Doppler absorption depth or Beer's length given in ^. Two separate mechanisms can change 
these familiar absorptive properties. These are self-induced transparency (SIT) which arises from 
dynamic nonlinearities associated with a strong and coherent pulse, and electromagnetically in- 
duced transparency (EIT) which arises when a second laser pulse interacts with the same excited 
state, thus inducing two-photon coherence between two different ground states. 

Potential physical realizations of such a system include the D line transitions for Rubidium or 
Cesium. For the D2 line of Rubidium, which we use as our model in this paper, the excited state 
lifetime is around 30 ns and the Doppler lifetime for room temperature vapor is about 0.5 ns, while 
the Beer's length is approximately 1 cm. We assume that the laser pulses are linearly polarized 
with a bandwidth greater than the hyperfine splittings of the excited state, which is approximately 
500 MHz, but with bandwidth narrow enough to resolve the two hyperfine ground states which are 
split by approximately 10 GHz. This implies the effective dipole moments, di and d2, are equal 
which means that taking fia = fJ-b is quite accurate. Our calculations use pulses with temporal 
duration around 2.5 ns, which is consistent with all of our approximations and assumptions. 

III. BRIGHT - DARK STATES 

The three-level MB model permits a useful parameterization in terms of bright and dark states 
[l7[, which help explain the interference effects caused by ground state coherence. We define the 
Bright and Dark states as: 

\B)^-^{na\l) + nb\2)) and |D) ^ J- ([7,*|1) - J]:|2)) , (5) 

where Qt = (l^^aP + l^bP)"*^'^ is the "total" Rabi frequency. In terms of the bright and dark states 
the Hamiltonian in Eq. ([1]) can be written simply as: 

H = hA\3){3\-^\B){3\-^\3){B\. (6) 

The interaction terms in H clearly depend only on \B). The orthogonal dark state \D) does not 
participate in the temporal dynamics so once population enters the dark state it becomes trapped. 



unless it evolves as a result of propagation. As clearly seen in the definition, the bright and dark 
state basis is a fully coherent, pure-state superposition of the two ground states. This basis helps 
provide a clear understanding of phase-coherent effects on two-pulse propagation. 

To more clearly understand these effects it is useful to convert the MB model into the bright- 
dark basis. This was done by Fleischhauer and Manka 123], and we will follow their lead. We 
will take a pure-state approach in this particular section, to focus on purely phase-coherent effects. 
Assuming a pure-state wavefunction of the form \ip) = ci|l)+C2|2)+C3|3), the probability amplitude 
equations in the original atomic basis are: 

ci = z— C3 (7a) 

C2=l — C2, (7b) 

C3 = i^ci+i^C2-iAc3, (7c) 

where the dot refers to d/dT. One can recover the density matrix given in Eqs. ([2]) by taking 
p = \'ip){'ip\. From Eq. ([5]) we can calculate the bright and dark probability amplitudes in terms of 
ci and C2- They are: 

CB = -^{n:ci + nic2) (8a) 

CD = — (rifeCi - flaC2)- (8b) 

For simplicity we will assume that the fields are unchirped and real, giving f^a = J7* and Qb = ^l- 
This allows us to rewrite the Eqs. ([7]) as: 

cb = "^^CD + «^C3 (9a) 

CD = -i-^CB (9b) 

C3 = i — CB - iAc3, (9c) 

where 

^D = ^{^A - ^fe^a) (10) 

is called the the dark Rabi frequency. We can write the corresponding Maxwell's equations for the 
total and dark Rabi frequencies as: 

-Q^ = -if^icBC^) (11a) 

d^D „ d /CDCI\ 



where //„ = /^fe = /^ as previously noted. 

A surprising result occurs if the pulses are matched temporally, i.e. (la/^a = ^b/^bi which 
gives VL£) = 0. This allows us to write Eqs. ([9]) and ([TT]) as: 

CB = ^-^C3 (12a) 

C3 = i—CB - iAc3 (12b) 

and 

-Q^ = -i^CBcl). (13) 

We see that all dark-state effects disappear, and Eqs. (J12p and (jlSh are identical in form to the 
Maxwell-Bloch equations for a single "total" pulse interacting with a two-level atom. Thus if the 
dark state is initially unpopulated and the pulses are matched temporally, then the population is 
constrained to states \B) and |3), and the three- level behavior is identical to that of a two- level 
atom. This explains the three-level simulton solutions derivation [14]. However we will find that 
small fluctuations prohibit such isolation of the dark state over long propagation distances. 

Kozlov and Eberly (KE) also used this fact to examine SIT-type vs. EIT-type transparency 
[20|. In SIT it is nonlinearities that cause the medium to be transparent even while the dark 
state is unpopulated and strong interaction occurs between the pulse and the medium. Thus SIT 
propagation occurs when |c_Dp ~ 0. As just shown, this implies that the three-level equations can 
be reduce to two-level form, and transparency only occurs if the two pulses are matched. In EIT it 
is the presence of a second pulse and population trapping in the dark state that cancels absorption 
of the pulses by the medium. Thus KE defined EIT propagation to occur when |c£)p ~ 1. In this 
case, the medium is transparent because the combined pulse-medium system is in the dark state, 
cancelling pulse-medium interaction. We will use this same definition throughout this paper. 

IV. MIXONIUM ANALYTICAL SOLUTIONS 

We wish to analyze two-pulse propagation through a partially coherent A system, i.e. a medium 
prepared with pij ^ but \pij\'^ < pi 1/0224°^ i ¥^ J- We do this by solving Eqs. ([2]) and ([3]) by using 
the Park-Shin (PS) Backlund method 28|]. We take the initial density matrix of each atom to be 
in a partially-coherent superposition of the two ground states, which we write in explicit form, for 



real a and /?, as 
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(14) 



where a^ and /3^ are the populations of ground states 1 and 2 respectively (we will always take the 
case a^ > 0^), and (f) is the phase of the partial coherence. We introduce a coherence parameter A 
that takes values between and 1, between a completely mixed state and a pure state respectively. 
We will assume the fields are tuned to the center of the inhomogeneous line for the two transitions. 
To find the solution, first we diagonlize p^^' in Eq. (J14p with the rotation matrix 



S 



I 



\ 



cost' 

sin^e" 





oi<t> n\ 



smt/e"'^ 
cos 6* 
1 



(15) 



where 



cos 9 



sm( 



C-P' 



V(C - /32)2 + A2a2/32 

—\a(3 
V(C-/32)2 + A2a2/32' 



(16a) 
(16b) 



and 



C = ^{l + [l-4(l-AV/3¥/n, 



(17) 



where C, and 1-^ are the positive non-zero eigenvalues of the matrix p^'^' . The rotation matrix S 
has an additional degree of freedom in that the first and second columns can be multiplied by an 
arbitrary phase, while still diagonalizing p^^'. In some plots of our analytic solutions we will make 
use of this fact. 

One can verify that S diagonalizes p^^' in Eq. ()14p through the operation 



(18) 





1-c 

\o oy 

We can solve the much simpler problem of pulses propagating through a medium prepared 
with initial density matrix given by Eq. (jlSp by using the Park-Shin J28l] Backlund transformation 



technique, just as we did in a previous paper [23|. These pulse solutions are: 

-1 



T 



n 



(d) 



2cosh( Ckz] +exp [ - + kZ{3C-2] 



2 cosh ( kZ{1 - C) ) + exp ( - + kZ{1 - 3C) 



(19a) 
(19b) 



where r is the nominal pulse width, the superscript (d) is to remind us that these are the solutions 
for two pulses propagating through a medium with a diagonal density matrix, and k/c is the inverse 
absorption length given by: 

2 y_» 1 + (Ar)2 ^ ' 

In the limit where Tg <^ r, the scaled inverse absorption depth 2k/c becomes the inverse Doppler 
absorption depth such that 2k/ c -^ an, which we previously defined as an = \/vr/2/iT2*/c. 



The solutions to the density matrix elements for a particular value of the detuning A are given 



by: 






1 



1 + (Ar)2 

(d) _ 1 

^22 - 1 + (Ar)2 

1 



C[|/iiP + (Ar)2] + (l-C)|/i2p 
CI/i2p + (l-C)[|/22p + (Ar)2] 



(d) 
P33 



l + (Ar 

(d) _ 1 

^12 - 1 + (Ar)2 

(d) _ 1 

^13 - 1 + (Ar)2 

1 



2, Cl/l3|' + (1-01/231' 



(d) 
P23 



1 + (Ar) 



C(/ii - iAr)/i2 + (1 - C)(/22 + ^At)/i2 
C(/ll - iAr)/i3 + (1 - C)/l2/23 
C/r2/l3 + (l-C)(/22-iAr)/23 



where the functions fij are: 

/ii = j 2 sinh(r/r - C^Z) - exp [T/t + (3C - 2)kZ] \ / D{Z, T) 

/22 = j - 2 cosh [T/t - CkZ) + exp [T/t + (3C - 2)kZ] \ / D{Z, T) 

f^2 = 2e^/^-^^-^^''^/D{Z,T) 
hs = 2i/D{Z,T) 
/23=2ze(2C-i)«^/^(^^r), 

and the denominator function D{Z, T) is given by: 

L>(Z, T) = 2 cosh (T/r - QkZ) + exp [T/t + (3C - 2)kZ] . 



(21a) 
(21b) 
(21c) 
(21d) 
(21e) 
(21f) 

(22a) 
(22b) 

(22c) 
(22d) 
(22e) 

(23a) 



Now that we have the solutions for a medium initially in the diagonal state, we obtain the 
mixonium solutions through the operations 
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(25) 



Eqns. (j24p and (j25p are the exact solutions to Eqns. ([2]) and ([3]) for a medium initially prepared 
in a mixed-state coherent superposition of the two ground states, as given in Eq. (J14p . One can 
see that the pulse and density matrix solutions are simply a rotation of the completely mixed-state 
solutions that we have previously solved 



23], with the eigenvalue of the initial density matrix 
taking the place of the ground state population. The invariance of the MB equations under such 
operations, which allows us to obtain the mixonium solutions, is discussed in further detail by Park 



and Shin 



28l |. It also complements previous numerical work which demonstrated the applicability 



of such "dressed- field" pulses [l6[] . 



PHASEONIUM ANALYTICAL SOLUTION ANALYSIS 



The pulse and density matrix solutions presented are clearly quite complicated. However, we 
can begin to understand them by starting with the pure state case where A = 1. In that case Eq. 
()17p becomes C = 1 ^iid Eq. (J16p simplifies to cos 6 = a and sin 9 = —f3. The pulse solutions then 
simplify to: 



^a = aOi<^) - Pe^'t'nf^ 



VLk = f5e-''^n^^'^ + a^ 



(d) 



(26a) 
(26b) 



The atomic solutions also simplify substantially. Since we are considering pure states, we can 
consider just the wavefunction. Thus Eq. (j25p can be simplified to: 

IV) = 51^^^)= pe-'^ a cf , 
V 1/ V^J 



(27) 
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where the diagonal probability amplitude solutions are given by: 



cT' = [fu - iAr] (28a) 

' Vl + (Ar)2^' ^ ^ ' 

-^^^ ^ ifi2 (28b) 



.(d) _ 1 



which can be derived from Eqs. (j21ap - (j21cp up to an overall phase. Thus the probability 
amplitudes on line-center (taking A = 0) can be written in a substantially simplified form as: 

ci = a/ii - /3/i2 (29a) 

C2 = (3fu + a/i2 (29b) 

C3 = /i*3 (29c) 

where we have assumed the phase of the ground state coherence = for simplicity, and we will 
do so for the remainder of this paper. Eqns. (|26p and (|29p mark a substantial simplification of the 
general mixed state solution. 

A. Phaseonium Input Regime: —kZ ^ 1 

We can further simplify the analysis by dividing the evolution into three distinct regimes. We 

will consider the asymptotic "input" as regime I, the asymptotic "output" as regime III, and the 

transfer zone in between as regime II. We study regime I by taking —kZ ^ 1. In this limit the 

pulse solutions become: 

2 fT \ 

Qa — *■ ct— sechl kZ (30a) 

T \t J 

Qb^P- sechl kz\ (30b) 

and the line-center probability amplitudes are given by 

ci — > atanh I kZ j (31a) 

C2^ P tanh ( kZj (31b) 

C3 — > —i sechl kZ | . (31c) 

These pulse and amplitude solutions are exactly the matched simulton solutions of Konopnicki and 
Eberly [ij] moving with group velocity Vg/c = (1 -|- kt)~^. This is an example of SIT type propa- 
gation, because |c£)p = in the limit, and the excited state is fully populated at the pulse peak. 

11 



However because of the particular pulse shape and their matching, the medium is transparent, but 
with a slowed velocity, possibly much slower than c. 

B. Phaseonium Output Regime: kZ ^ 1 

Similarly the output regime III can be considered by taking kZ ^1. In this limit the pulse 
solutions are: 

Qa^ -/?- sechf- J (32a) 

fife^a- sechf- j, (32b) 

and the line-center probability amplitudes are simply the constant values 

c\ — > —a (33a) 

C2 ^ -/3 (33b) 

cg -^ 0. (33c) 

Just as in the input regime I, these pulse solutions are also matched simultons but now moving with 
the vacuum velocity c. One can see from Eqns. (f33l) that the excited state probability amplitude 
is 0. This is because all population is now in the dark state, such that \cj:)^ = 1 in the limit. 
Thus the ground state amplitudes become constant and EIT type propagation occurs, allowing the 
pulses to propagate without any interaction with the medium. We note from Eqs. (j32p that the 
pulse shapes are matched, and with ratio 



which is exactly the ratio predicted by the Dark Area Theorem [ic 

As remarked previously when discussing the Hamiltonian, we know that the dark state is com- 
pletely decoupled from the dynamics. Thus if during the propagation the pulses become matched 
to the medium, as in (j34p . no further population transfer can occur. We see that when A = 1, 
meaning the initial state of the density matrix is a pure state, the analytic solutions describe a 
transfer of SIT-type simulton pulses to EIT-type pulses in the dark state, where the population is 
trapped and remains constant. We also see from these analytic solutions, that the original simul- 
ton solutions [IJ] are simply limiting cases of our more general solutions. The simulton solutions 
propagate in stable form without any change and without ever becoming trapped in the dark state. 
However this picture is not complete as we now see that transfer to the dark state still occurs. 

12 



We plot the analytic pulse solutions, JIq and Jl;,, given in Eqs. ([26]) . in the left frame of Fig. 
[21 and the line-center excited state probability, |c3p, given in Eq. (j29p . in the right frame of Fig. 
[2j The numbers in each frame correspond to a particular time point (note that the time points 
are chosen to illustrate interesting changes that are occurring and are not uniformly spaced). 
Examining both of these figures together we clearly see that in the initial simulton regime I, the 
pulses propagate in the bright state, causing excitation into the excited state (frames 1 and 2). Then 
in the transfer regime II (frames 3-5), we see the relative magnitude of the pulses changing, and 
the phase of the Stokes pulse changing sign, along with a decrease in the excited state probability. 
Finally, in regime III (frame 6) we see both pulses propagating without change and the upper level 
is not excited at all, since the pulses and medium are now in the dark state. 
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FIG. 2: Plots of the analytic pulse solutions given in Eq. ([SG]) on the left, and of the analytic excited state 
population solutions given in Eq. (|29p on the right. The pump pulse, U,a is the solid curve, and the Stokes 
pulse, ilf, is the dashed curve. The horizontal axis is x in units of k/c. The vertical axis is the pulse Rabi 
frequency in units of r~^ (left) and the excited state probability P33 — |c3p (right). The background is 
slightly shaded to indicate the presence of the lambda medium. The solid curve is the pump pulse, U,a, and 
the dashed curve is the Stokes pulse, ilf,. The plot shows the simulton transfer process as an exchange process 
between bright-state simultons to dark-state simultons. The excited state is heavily populated during the 
input propagation regime, when the pulses and atoms are in the bright state. However as the pulses and 
atoms transfer into the dark state, the atoms no longer absorb any pulse energy, and the medium becomes 
transparent. Parameters: a^ — 0.8, /3^ — 0.2, r = STj*; £^nd A = 1.0. 
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VI. MIXED STATE ANALYSIS 

We now consider the pulse and density matrix solutions for a medium prepared in an arbitrary 
mixed initial state (i.e., < A < 1). From Eq. (j24p . we obtain the individual pulse solutions which 
are given by: 

Qa = cos 0Oi^) + sin en'^^ (35a) 

Qfe = - sin 0Qi<^) + cos en'^^^ , (35b) 

where cos 9 and sin 6 are given in Eq. p6|) and ^a and Qj^ are defined in Eq. (|19|) . The mixed- 
state density matrix solutions are given in Eq. (j25p and are clearly quite cumbersome for this 
general mixed state case. The simplest of these, and the one which provides the most insight into 
the physics is the excited state probability, p^^. Written explicitly the line-center solution is: 

P33 = : ^ : -T2- (36) 

[2 cosh(r/r - CkZ) + exp {T/t + (3C - 2)kZ)] 

Once again we will examine these solutions in the input and output regimes to help understand 

their underlying features. 



A. Mixonium Input Regime: —kZ ^ 1 

In the input regime I, by taking the limit —kZ ^ 1, we find: 

2 /T \ 

Q,a -^ cos 6— sechi C,kZ J (37a) 

2 fT \ 

Qi) —^ — sin 6— sech I ^kZ J , (37b) 



while the excited state probability is: 



P33^Csech2(^-CA€Z). (38) 



Just as in the pure state case, Eqns. (I37p are matched simulton pulses. However they are now 
generalized to non-pure states. The most immediate differences between the two cases are the 
modification to the pulse amplitudes, and that the excited state no longer reaches a value of 
P33 = 1 at the pulse peak, but rather P33 = (". This results in a group velocity of the pulses, 
Vg/c = (1 + (^Kr)~^, which is diff'erent by the presence of the factor C relative to the pure state 
case. This indicates that the interaction between the pulses and the medium is directly affected by 
the medium's mixed-state nature, and that the interaction strength is governed by the parameter 
C, which is simply the eigenvalue of the initial mixonium density matrix. 
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B. Mixonium Output Regime: nZ ^ 1 

We look at the mixed-state output regime III by taking kZ ^ 1. In this hmit the pulse solutions 
are: 

Ha -^ sin 6*- sech( (1 - C)kZ j (39a) 

0(,^cos6'- sechf {I - C)kz\ (39b) 

and the line-center excited state probability is now: 

P33 ^ (1 - C) sech^ (^ - (1 - 0^^Z\ . (40) 

The output pulses are again matched simultons, and have quite similar features to the output pulses 
for the pure state. The major difference is that the dark state can no longer be fully populated so 
interaction between the pulses and medium continues, however because of SIT the pulses continue 
to propagate without absorption. The group velocity is not c, rather it is Vg/c = [1 -|- (1 — C)i^t]~^ , 
and at the peak of the pulse the excited state probability is P33 = 1 — C instead of 0. 

Both the input and output regimes show modified excited state populations when compared 
to the pure state case, and the modification in both regimes is governed by the same parameter 
C. We will refer to this parameter (" as the interaction parameter. Its value gives the maximum 
population of the dark state, and thus determines how strongly EIT can cancel the pulse-medium 
interaction. In the pure-state case C = 1, and all interaction ceased once the pulses were in the 
dark state. However as just shown, when the medium is initially in a mixed-state, we have C < 1) 
and interaction continues. The value of (" decreases until the limiting completely mixed state is 
reached, where A = 0, giving (^ = a^ (for a^ > /3^). In Fig. [3] we plot the value of this interaction 
parameter as a function of the coherence parameter A for a variety of initial medium preparations. 
The plot shows the range of the parameter (^ and its dependence on initial medium preparation. 
However as the medium approaches the pure-state case, all curves converge, and thus all interaction 
is cancelled no matter how the population is initially distributed. 

Next we plot the pulse solutions for a mixed-state medium with parameters a^ — /3^ = 0.6 in 
Fig. m The left figure corresponds to A = 0.8 and the right figure corresponds to A = 0.2. The 
solutions for the excited state population are plotted underneath in Fig. [5l For the pulses we 
see a very similar propagation behavior to that of the pure-state case plotted in Fig. [2j Aside 
from a slight difference between the relative pulse amplitudes the plots look very alike. The main 
difference is seen in the plots of the excited state population. The input pulses shown in frame 1 
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FIG. 3: Plots of the interaction parameter defined in Eq. ([TT]). The horizontal axis is the coherence 
parameter A, and the vertical axis is the interaction parameter Q. We plot the interaction parameter for a 
variety of media preparations ranging from a"^ — 0^ ~ 0.2 up to 0.8. It ranges in value from ( = a^ for A = 
(assuming a^ > /9^) to C = 1 for A = 1. 

no longer cause complete excitation into the excited state so ^33 7^ 1 even at the pulse peak. More 
importantly the output pulses are no longer completely decoupled from the medium, since P33 ^ 
in the output regime III as seen in frames 5 and 6 of Fig. [5l This feature is more enhanced as the 
value of A gets smaller. The fact that the medium is not in a fully coherent pure-state superposition 
implies that the dark state cannot be fully populated and thus the pulse and atoms still interact. 
However because of SIT the medium is still transparent to the pulses. 

VII. PULSE AREA 



We showed in section IIIII that in a pure-state medium, with the pulses initially matched, such 
that the dark state is cancelled, the medium will behave exactly as a two-level medium when 
described in the dressed bright-dark basis. In this basis the two pulses must be combined and 
thought of as one "total" pulse, with total Rabi frequency defined as Qt = (l^aP + l^feP)"*^ • We 
also found in the previous section that even though interaction occurred between the pulses and 
atoms in a mixonium medium, the medium still appeared transparent due to SIT. Thus we expect 
some elements of the two-level area theorem to hold for the bright pulse area. With the area of a 
pulse defined to be 



AiZ) 



n{Z,T)dT, 



(41) 
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FIG. 4: Plots of the analytic pulse solutions given in Eq. ([35]) for a mixonium niedium for two different 
values of A. The horizontal axis is x in units of k/c, and the vertical axis is the pulse Rabi frequency in units 



of 



The pump pulse, Qa, is the solid curve, and the Stokes pulse, fif,, is the dashed curve. The plot 



shows the siniulton transfer process as an exchange process between bright-state siniultons to dark-state 
simultons. Parameters: a^ = 0.8, f3^ = 0.2, r « 3T^, and A = 0.8 (left figure) A = 0.2 (right figure). 
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FIG. 5: Plots of the analytic excited state population solutions given in Eq. (|36p for a mixonium niedium. 
The horizontal axis is x in units of k/c, and the vertical axis is the excited state probability ^33 . The plot 
shows the excited state being heavily populated during the input propagation regime, when the pulses and 
atoms are in the bright state. However unlike the pure-state case, the excited state never reaches ^33 = 1. 
As the pulses and atoms transfer into the dark state the pulses and atoms have a weaker interaction, however 
the interaction is never completely cancelled and P33 7^ 0. Parameters: a^ — 0.8, /3^ — 0.2, t « B^j*, and 
A = 0.8 (left figure) A = 0.2 (right figure). 
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the pulse areas of the general mixed state solutions (j35p can be shown to be 



, , „, / cos 6 sin ^ \ , , 

A rrys ^ (-Sm.0 COSO \ /-^i N 



where 



/i(Z) = v/l + e2(2C-i)«^. (43) 

A remarkable results occurs when considering the total-pulse area: 

/oo 
-oo 

Since the analytical solutions are temporally matched (i.e. d/dT{Vta/^b) = 0), the total pulse area 
can be written as the sum of the squares of the individual pulses areas. Thus the area of the total 
Rabi frequency from the analytic solutions is simply: 



At{Z) = ^AI{Z) + AI{Z) = 2^. (45) 

This indicates a remarkably broad connection to SIT in which the pulse area must obey the area 
theorem: 

Even for an arbitrary mixed-state medium, the total pulse area remains constant and equal to 27r. 
We will elaborate on this SIT connection in the next section via numerical solutions. 

VIII. NUMERICAL SOLUTIONS 

We now focus our attention on examining the consequences of both the modified dark state 
properties of a mixed state medium, and the connection of our pulse solutions to SIT. We do this 
by testing our exact analytical solutions in a more realistic experimental setting, through numerical 
solutions to the MB equations. This allows us to test the general utility of insights suggested by the 
highly specialized exact matched sech shaped analytical input pulse solutions for infinite medium 
length. We will use gaussian input pulses defined as: 

0(-) = ^e-S and l^f ^ = ^e'S, (47) 
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and "super-gaussian" input pulses defined as: 

0M = ^e-(^ and f^f ^ = ^e'c^, (48) 

where Aa gives the pump pulse area as defined in ()4ip . Tq is the nominal pump pulse width, 
r(l/4) ~ 3.6 is the gamma function, and similarly for the Stokes pulse. We replace the infinite 
uniform medium by a medium with definite entry and exit faces. 

We will first review pulse propagation in phaseonium, a pure-state medium. We will show how 
the SIT-like propagation predicted in Sec. IIIII can be realized by using matched input pulses with 
ratio ila/^b = o//3 such that |cdP = 0. However small fluctuations eventually cause these pulses 
to reach the dark state. Next we will show numerical solutions with unmatched input pulses that 
demonstrate pulse matching, similar to results shown in ref. [20]. We will contrast these results 
with solutions for pulse propagation in a mixonium medium. Because the dark state can never be 
fully populated for pulses propagating in mixonium, SIT always plays a role, and EIT dominance 
is weakened. The difference is most pronounced for long media that are many absorption depths 
long, however even short media show some of the characteristics. 

A. Phaseonium - Matched Input Pulses 

As shown in Sec. IIIH when the pulses are initially matched, with ratios Qa/^b = ck//5, such that 
|cDp = 0, the three- level A system reduces to that of a two- level system regardless of the composite 
bright-pulse shape. We plot such a solution in Fig. [U where we use matched super-gaussian input 
pulses as defined in Eq. ([18|) . with equal pulse widths Ta = t^ = STg* and areas Aa = 2.3Q7r 
and Ah = 2.3/?7r. The medium is prepared in a pure state such that A = 1.0 with populations 
a^ — 0^ = 0.6. This solution clearly illustrates the three regime language that we have introduced. 
In regime I (frames 1-3) the pulses act as simultons, and the pulse-medium interaction is exactly 
that for resonant two-level pulse propagation when considered in the dressed state basis. Then 
in regime II (frames 4-5) the pulse amplitudes rapidly change leading to the dark-state output 
regime III (frame 6). The solution exhibits SIT type propagation in regime I, where \c]:)\^ ~ 0, but 
beginning near frame 4 the pulse amplitudes begin to change, and eventually end up with ratio 
ilfl/ilfe = — /3/a in frame 6. Here all population is in the dark-state and |cdP ~ 1, and EIT type 
propagation takes over. 

The initial SIT type behavior is further illustrated in Fig. [7] where we plot the areas of the 
individual pulses, as well as the area of the bright pulse. After a few absorption depths, the total 
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FIG. 6: Plots of numerical pulse solutions of Eqs. ^ for a pure-state phaseonium medium, with super- 
gaussian input pulses. The horizontal axis is x in units of k/c, and the vertical axis is the pulse Rabi 
frequency in units of r^^. The solid curve is the pump pulse, fJa, and the dashed curve is the Stokes pulse, 
57b- The plot shows matched input pulses with ratio fla/flb — ck//3, such that |cdP ~ reshaped just as in 
normal SIT. However after some propagation distance their ratios change to fla/^b = ^f^/a as predicted by 
the dark area theorem so that |cdP ~ 1. Parameters: a^ = 0.8, /3^ = 0.2, Ta = t^ ~ t Ri 3r2*j ^a = 2.3a7r, 
Af, — 2.3/37r, and A = 1.0. The value of ^ is chosen so that Vg/c — 1/2 inside the medium. 
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FIG. 7: Numerically integrated area of the individual pulse areas as well as the total Rabi frequency area 
for the pulse solutions shown in Fig. [B] The vertical axis is the pulse area. The solid curve is the area of 
the pump pulse, the dashed curve is the area of the Stokes pulse, and the dot-dashed curve is the area of 
the total Rabi frequency. The medium initially behaves as a two-level medium, and the bright area changes 
to 27r. However after some propagation distance the pulses transfer to the dark state, while the bright area 
remains constant through this change. 
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pulse area is quickly changed from its input area of At = 2.3-k to 27r. Then after propagating as 
simultons for a while, the pulses quickly change to the dark state. While the individual pulse areas 
are rapidly changing, we see completely static behavior for the bright area, confirming our analytic 
result in Eq. (031). 

After the initial reshaping caused by the two-level SIT behavior, the pulses propagate for some 
time as simultons. The simulton solutions are exact solutions, so the reason for the pulses to 
transfer to the dark-state does not readily present itself. Only perfect matched sech shaped pulses 
remain as simultons. Any small perturbation from these exact shapes implies |c£)p > 0. While the 
initial dark-state probability may be small, any non-zero value will eventually lead to the dark state 
if the propagation distance is long enough. In this particular example, the small perturbations are 
caused by numerical roundoff error. The analytic solutions actually predict this behavior. It is the 
exponential term in the denominator of Eqs. (jlOp which describes deviations from the simulton 
shape. Thus we see that simulton pulse propagation is in fact an inherently unstable propagation 
scenario, and that small fluctuations will always lead to the dark state. 

B. Phaseonium - Mismatched Input Pulses 

Next we examine what happens if the pulses are initially temporally mismatched, but where 
the medium is still prepared in a pure state. We show the numerical pulse solutions in Fig. [HI 
where we take the medium to be prepared with a^ — /3^ = 0.6, and the pulses to be gaussian shape 
with width r^ = Tj,/2 = STg , a temporal mismatch of 2T(j, and areas Aa = 1.27r and A^ = O.Stt. 
In this case, the pulses do not propagate as simulton pulses, since they are not initially matched 
(i.e. Qa/^b 7^ constant). However as they propagate they are quickly reshaped into two-peaked 
but matched pulses with ratio i^a/^b = — /5/ct, such that the dark-state is fully populated and 
BIT type propagation occurs. This re-shaping occurs because the dark Rabi frequency is nonzero 
for mismatched pulses. In this example the medium is only 10 absorption depths long, and most 
re-shaping is complete after about 5 absorption depths. This plot confirms previous KE results 
showing very similar behavior [20t]. 

We also plot numerically integrated pulse areas in Fig. [9] for both the individual pulse areas as 
well as the area of the total Rabi frequency. This plot further illustrates the distinction between 
the matched and mismatched cases. The pulse areas are quickly modified as the pulses match, but 
as soon as they reach the dark state, the pulse areas become constant. The SIT area theorem does 
not apply in this case, and no prediction is possible as to the final pulse areas, unlike the previous 
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FIG. 8: Plots of numerical pulse solutions of Eqs. ([2|) for a pure-state phaseoniuni medium, with gaussian 
input pulses. The horizontal axis is x in units of k/c, and the vertical axis is the pulse Rabi frequency in units 
of r^^. The solid curve is the pump pulse, fi^, and the dashed curve is the Stokes pulse, f2(,. The plot shows 
mis-matched input pulses quickly reshaped into matched pulses with ratios Qa/^b ~ ~P/oi as predicted by 



the dark area theorem and the analytic solutions. Parameters: a^ — 0.8, 0^ — 0.2, Tq — Tb/2 — t k, 3Tr 
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FIG. 9: Numerically integrated area of the individual pulse areas as well as the total area for the pulse 
solutions shown in Fig. [S] The horizontal axis is x in units of k/c, and the vertical axis is the pulse area. 
The solid curve is the area of the pump pulse, the dashed curve is the area of the Stokes pulse, and the 
dot-dashed curve is the area of the total Rabi frequency. The pulses initially reshape as they match, however 
after a short propagation distance the pulses enter the dark state, and the pulse areas become constant. 
The SIT area theorem does not apply in this case. 
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example. 

These two examples of matched and mis-matched input pulses highlight a feature common to 
pulse propagation in phaseonium. That is, the dark-state always dominates. Pulses will always 
end up matched and always with a ratio that satisfies the dark area theorem, so that absorption 
is cancelled and EIT plays a dominate role. However, for matched input pulses with no dark state 
population, two-level physics initially describes the propagation, and the dark-state dominance 
takes much longer to appear. While we can identify the input, output and transfer regimes in the 
matched input example, we cannot do the same for the mismatched input pulses. We will see in the 
next section that the mixonium medium modifies the absorptive properties causing the dark-state 
dominance to be replaced with SIT like effects. 

C. Mixonium - Mismatched Input Pulses 

We now examine the effects that mixonium has on mismatched pulse propagation. Matched 
input pulses (no matter the shape) with input ratios given by the analytic solutions in Eqs. (j37p 
will behave in a similar manner to the matched pulse example in the previous section, so we do not 
plot the results here. The difference is simply that the dark state can no longer be fully populated 
as discussed in Sec. Wl\ and thus the pulses continue to cause excitation into the excited state. 

The inability of the dark state to be fully populated in mixonium has a profound effect on 
mis-matched pulses propagating through many absorption depths (the result is less dramatic for 
short media). We plot the pulse solutions for the same parameters as the solutions plotted in the 
previous section, with a^ — /3^ = 0.6, gaussian pulse shapes with duration r^ = rf,/2 = 3T2* and 
an offset of 2ra, and pulse areas of Aa = 1.2tt and At, = O.Svr, except we take the medium to be in 
a mixed state with A = 0.8. We plot these pulse solutions in the left frame of Fig. [lOl We see a 
very similar behavior to the previous solutions with one almost unnoticeable difference. That is, 
the output pulse ratio is given by ila/^b = tan ^ = — -^"/5/(C ~ /3^) as predicted by the mixed-state 
analytic solutions. However, when we plot the areas of the pulses in the left frame of Fig. [11] we 
notice now that the total pulse area is very close to 27r, in contrast to the pure-state solution shown 
in Fig. m In fact in this example, the pulses were still being reshaped, thus we will examine what 
happens if the medium is slightly longer. 

The difference between phaseonium and mixonium propagation is greatly magnified when the 
pulses propagate through a longer medium. We plot these solutions in the right frame of Fig. [10] 
which show continued reshaping of the pulses until they are matched sech shaped pulses, exactly as 
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FIG. 10: Plots of numerical pulse solutions of Eqs. ^ for a mixed-state medium, with gaussian input 
pulses. The horizontal axis is x in units of k/c, and the vertical axis is the pulse Rabi frequency in units 



of ■ 



The left and right frames are identical except for different medium lengths. The solid curve is the 



pump pulse, fla, and the dashed curve is the Stokes pulse, fit- The plot shows mis- matched input pulses 
quickly reshaped into matched pulses, similar to the pure state case. Parameters: a^ = 0.8, /3^ = 0.2, 
Ta = T(,/2 = r w 3T2*, temporal mismatch of 2ra, Aa — 1.2tt, At, — O.Stt, and A — 0.8. 
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FIG. 11: Numerically integrated area of the individual pulse areas as well as the total Rabi frequency area 
for the pulse solutions shown in Fig. 1101 The horizontal axis is x in units of k/c, and the vertical axis is the 
pulse area. The solid curve is the area of the pump pulse, the dashed curve is the area of the Stokes pulse, 
and the dot-dashed curve is the area of the total Rabi frequency. We see initially rapid change in the pulse 
areas as they are matched, consistent with pure-state behavior. However because the dark state can never 
be fully populated the medium always behaves as an SIT like medium, and eventually the bright pulse area 
changes to 2tt area. 
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FIG. 12: Plots of numerical excited state population solutions of Eqs. ([2|) for a mixed-state medium. Each 
frame corresponds exactly to the same frame on the right hand side of Fig. 1101 The horizontal axis is x in 
units of k/c, and the vertical axis is the excited state population. The plot shows that even after the pulses 
are matched, the dark state can never be fully populated in a mixed-state medium, and thus the excited 
state continues to be populated. Parameters: a^ = 0.8, 0^ = 0.2, Ta = Tb/2 = t k, STj*, temporal mismatch 
of 2Ta, Aa = 1.27r, Afc = O.Svr, and A = 0.8. 

given by the output analytic solutions in Eqs. (j39p . We plot the pulse areas for this example in the 
right frame of Fig. [11] and we see that after the initial rapid reshaping which matches the pulses, 
the areas of the pulses continue to change until the bright pulse area reaches Ai, = 27r. The modified 
interaction properties of mixonium cause SIT to dominate and reshape the pulses. Thus the EIT 
dominance that was exhibited in the pure-state case, is now replaced with SIT dominance for 
mixed-state media. For this same long medium example we also plot the excited state population 
in Fig. [T2j We see that even after the pulses are matched, the excited state can never be completely 
decoupled. Thus, the dark state can never be fully populated, with its maximum value given by 
the interaction parameter C,. The inability of the dark state to be fully populated allows SIT to 
continuously reshape the pulses until they agree with the analytic solutions. 

IX. CONCLUSIONS 

We have presented new solutions of the Maxwell-Bloch equations, both analytic and numerical, 
applicable to a "mixonium" medium, where the term mixonium implies a A medium prepared in a 
partially phase-coherent superposition of the ground states. This medium offers a new contrast to a 
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pure "phaseonium" medium, where the ground states are prepared in a completely phase-coherent 
superposition of the two ground states. The partially coherent medium is experimentally realistic, 
whereas pure-state preparations are difficult to achieve. 

The analytic solutions for the pulses and density matrix elements were obtained via the Park- 
Shin Backlund transformation method 
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28( 1 . Consistent with our previous work 
identified three distinct regimes of interest for the analytic solutions. For the pure-state case, we 
identified our solutions in the asymptotic input regime to be equivalent to the well known simulton 
solutions [14^ . Our analytic solutions then describe the transfer of these input simulton pulses to 
simulton solutions completely in the dark state in the output regime. We identify this behavior 
as the transfer of pulses propagating with completely SIT like behavior to completely EIT like 
behavior, where we use the same definition of SIT and EIT like behavior as given in Ref. [201]. 
Using numerical solutions to the Maxwell-Bloch equations we were able to show this SIT simulton 
to EIT simulton behavior for pulses with matched but different shapes from the analytic solutions. 

The analytic solutions in the general mixed-state case, allow us to identify an interaction pa- 
rameter that determines the maximum population of the dark state. In the pure-state case, the 
dark state can become fully populated and all interaction between the pulses and medium is elim- 
inated. However in the mixed-state case the dark state can never be fully populated and thus the 
excited state cannot be decoupled. We studied the effects that this modified interaction parameter 
has on dark-state propagation dynamics by numerically solving the Maxwell-Bloch equations. In 
the pure-state case, our numerical results confirmed previous results that mismatched input pulses 
with gaussian shapes quickly match and propagate unchanged as the dark-state is fully populated 
prohibiting any coupling to the excited state [20]. In contrast, in the mixed-state case these same 
mismatched pulses never reach complete EIT type propagation since the dark state is never fully 
populated. Unlike the pure-state case, EIT type effects cannot completely cancel SIT type effects, 
and the pulses continue to be reshaped into matched sech shaped simultons, matching shape to 
our output analytic solutions. 

We have been able to demonstrate that the two-level McCall-Hahn area theorem still plays a role 
even in this three-level system. The composite or "total" two-pulse Rabi frequency of the analytic 
solutions has constant 27r area during all stages of propagation. SIT propagation effects cannot 
be completely cancelled for mixed-state propagation, causing the composite total Rabi frequency 
area to evolve toward 2tt just as in single pulse two-level SIT. Only in the pure-state case where 
EIT can completely cancel SIT, can this behavior be avoided. Thus we see that the two-level area 
theorem continues to play a role in two-pulse propagation through a three-level medium suggesting 
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a possible three-level area theorem for the total Rabi frequency. 
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